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Coherent controlization, i.e., coherent conditioning of arbitrary single- or multi-qubit operations 
on the state of one or more control qubits, is an important ingredient for the flexible implementation 
of many algorithms in quantum computation. This is of particular significance when certain sub¬ 
routines are changing over time or when they are frequently modified, such as in decision-making 
algorithms for learning agents. We propose a scheme to realize coherent controlization for any 
number of superconducting qubits coupled to a microwave resonator. For two and three qubits, we 
present an explicit construction that is of high relevance for quantum learning agents. We demon¬ 
strate the feasibility of our proposal, taking into account loss, dephasing, and the cavity self-Kerr 
effect. 
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I. INTRODUCTION 

The ability to coherently control and manipulate in¬ 
dividual quantum systems lies at the heart of modern 
quantum technologies and applications in quantum infor¬ 
mation [1-3] . Any quantum computation can be realized 
as a sequence of elementary quantum gates [4], which 
are highly-controlled quantum interactions of few qubits 
at a time, and quantum measurements. Prominent ap¬ 
plications include, e.g., quantum algorithms for efficient 
factoring [5, 6] and quantum simulation [7-11]. More 
recently, applications of quantum algorithms to certain 
problems in machine learning, including data classifica¬ 
tion [12] and search engine ranking [13, 14], have been 
proposed. Other recent proposals, which are of particular 
interest for the current paper, are the quantum-enhanced 
deliberation of learning agents in the context of quan¬ 
tum artificial intelligence [15, 16] and the notion of au¬ 
tonomous and adaptive devices for quantum information 
processing [17]. In parallel to these theoretical devel¬ 
opments, the design of experimental implementations of 
quantum computational architectures in systems such as, 
e.g., trapped ions [18-21] and optical setups [22-24] has 
been greatly advanced. 

In addition, the progress in controlling supercon¬ 
ducting (SC) quantum systems [25-28] has significantly 
strengthened the role of SC qubits (see, e.g., [29, 30]) as 
contenders for the realization of quantum computational 
devices. In particular, 1-D [31] and 3-D [32, 33] trans- 
mons — SC qubits impervious to charge noise — appear 
as promising candidates. Initial studies of complex algo¬ 
rithms and gate operations [34, 35], fault-tolerance [36- 
38], and hardware efficiency [39] in such systems further 
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raise the hopes for a scalable quantum computational 
architecture using SC circuitry. Moreover, coupling SC 
qubits via microwave resonators permits access to the do¬ 
main of cavity QED [40-42]. There, the exceeding level 
of control over the combined system can be utilized, e.g., 
to resolve photon number states [43], or to determin¬ 
istically encode quantum information in the resonator 
states [44, 45]. We shall draw from this rich quantum 
optics tool-box in the following. 

Notwithstanding these developments, a paradigm with 
respect to which all of the above implementations are typ¬ 
ically deemed successful is the realization of the unitaries 
of a universal set of quantum gates, see, e.g., Refs. [3, 46]. 
While a universal set of gates in principle allows any uni¬ 
tary to be efficiently approximated, there are some tasks 
for which this approach lacks a certain flexibility. For 
instance, when specific subroutines of an algorithm re¬ 
quire modifications in between individual runs. Promi¬ 
nent examples include the period-finding subroutine in 
Shor’s algorithm [6], which is typically used several times 
for different functions throughout the algorithm, or when 
the phases corresponding to different unitaries are to be 
estimated using Kitaev’s scheme [47]. This issue is also 
of particular importance when quantum subroutines are 
used in the decision-making of learning agents [48] , which 
update their subroutines based on experience gathered 
throughout the learning process. In all these cases, a 
set of unitaries, which are applied conditionally on the 
states of an ensemble of auxiliary qubits, is modified in 
subsequent applications of the subroutines. It is hence 
of significant interest to establish a method of coherent 
controlization — a mapping from a set of unitaries on a 
target Hilbert space to a single controlled operation on a 
larger (control & target) space — that is independent of 
the chosen set of unitaries. 

Unfortunately, it is not possible to design generic quan¬ 
tum circuits that achieve this conditioning independently 
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of the selected unitaries when only single uses of the uni- 
taries in question are permitted [49-51]. On the other 
hand, physical implementations of the unitaries that are 
to be controlled are typically already conditioned on fix¬ 
ing some degrees of freedom such as spatial locations (a 
laser beam illuminating an ion; an optical element being 
placed in the path of a light beam), or resonance frequen¬ 
cies. This practically allows the realization of purpose- 
built schemes that “add control” to unspecified unitaries, 
e.g., in optical setups [52], or trapped ions [48, 51]. 

Here, we propose a modular and adaptive implemen¬ 
tation of coherent controlization in a superconducting 
system of transmon [31] qubits coupled to a microwave 
resonator. In the dispersive limit, the coupling between 
these systems can be understood as well-resolved shifts 
of the cavity frequencies, dependent on the qubit states, 
or vice versa, shifts of the qubit frequencies conditioned 
on the cavity state. Based on this principle, our pro¬ 
tocol is assembled from unconditional displacements of 
the cavity mode and qubit operations conditioned on the 
vacuum state of the resonator, similar as in Refs. [44, 45]. 
We present a detailed construction of our protocol for two 
and three qubits, and we give a recipe for up-scaling our 
scheme to an arbitrary number of qubits. As a corner¬ 
stone of our investigation, we include an in-depth analysis 
of effects detrimental to the success of our protocol. For 
the strongest source of errors, the cavity-self Kerr, we 
provide analytical estimates of the disturbance, and dis¬ 
cuss methods to reduce it. Using numerical simulations 
to further take into account decoherence effects such as 
amplitude- and phase-damping, we hence show that our 
scheme for coherent controlization can be implemented 
using current superconducting technology. 

This paper is organized as follows. In Sec. II we briefly 
review the basic concepts for our proposal: A definition 
for coherent controlization is given in Sec. II A, illustrated 
by an application to learning agents in quantum artificial 
intelligence in Sec. IIB, before we give a short descrip¬ 
tion of the superconducting transmon qubits that we con¬ 
sider here in Sec. IIC. The conceptual centrepiece, the 
coherent controlization protocol for two superconducting 
qubits, is introduced in Sec. Ill, where we first discuss 
the idealized basic protocol, before turning our attention 
to the influence of the Kerr effect and decoherence. Fi¬ 
nally, we extend our protocol to three qubits and beyond 
in Sec. IV. 


II. FRAMEWORK 

A. Coherently controlling unknown unitaries 

At the heart of many quantum computational algo¬ 
rithms lie subroutines in which operations of choice on 
a finite-size register of qubits are performed condition¬ 
ally on the state of a control qubit. This is the case, for 
instance, in Kitaev’s phase estimation subroutine [47], 
where, the 2”“^-fold application of a phase ro¬ 


tation U{ip) is executed only if the n-th ancilla qubit is 
in the state |1)„. Similar subroutines feature also in 
Shor’s factoring algorithm [6]. When U{(p) is specified, 
the corresponding subroutine can be efficiently approx¬ 
imated by combining operations from a set of universal 
quantum gates (see, e.g., [46]). The decomposition into 
the universal gates, and their assembly to form the sub¬ 
routine clearly requires some (classical) computational 
effort along with (some) knowledge of U{(p). This be¬ 
comes a practical impediment when a device implement¬ 
ing said subroutine is to be used consecutively for differ¬ 
ent choices of U{ip), possibly diminishing any computa¬ 
tional speed-up with respect to purely classical devices. 
In particular, this issue is of crucial interest for the design 
of quantum-enhanced autonomous learning agents [15], 
where the quantum speed-up concerns the deliberation¬ 
time, and the agents need to update their subroutines 
throughout the learning process. 

It would hence be highly desirable to have access to 
fixed global operations, let us call them A and B, on the 
target and control registers which allow turning an un¬ 
specified local operation U{ip) into its controlled version 
Ctrl— U{(j)), such that AU{ip)B = Ctrl— U{(j)). However, 
this requirement cannot be met by any fixed A, B for all 
[/((/?), see [50], and therefore, in particular, generic sys¬ 
tem independent controlization is not possible when the 
action of U{ip) on the target Hilbert space is unknown. 
Fortunately, most practical realizations of (unitary) op¬ 
erations are not strictly local with respect to the target 
Hilbert space, but are already conditioned on some addi¬ 
tional degrees of freedom. For instance, for quantum in¬ 
formation encoded in photons, the optical elements must 
be placed in the beam path, conditioning the operations 
on spatial degrees of freedom. Another example are laser 
pulses driving transitions between qubit states encoded 
in ions, which must be at resonance, conditioning the 
transformations on the correct frequency. Such implicit 
conditioning on additional degrees of freedom can be ex¬ 
ploited to “add control” to unitaries that are unknown 
in the sense specified above [51, 52]. 

Following Ref. [48], we shall refer to mappings from a 
set of operations {Ui} on the target Hilbert space Ht to 
operations U on the joint Hilbert space He 0 Ht of the 
control and target systems, such that 

( 1 ) 

i 

for some (orthonormal) basis {|i)}c of He, as well as 
to any specific physical realization of such mappings as 
coherent controlization. In the following, we shall dis¬ 
cuss how general coherent controlization can be imple¬ 
mented in superconducting qubits, providing a modu¬ 
lar and adaptive architecture for quantum computational 
tasks. For a detailed analysis of the feasibility and scala¬ 
bility of our proposal, we provide a concrete example for 
an application of coherent controlization in the decision 
process of learning agents, where the adaptive character 
of our proposal is of particular relevance. 
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B. Coherent controlization in the context of 
learning agents 

In the model of projective simulation (PS) [16], an 
autonomous learning agent draws upon previous expe¬ 
rience to simulate its future situation in a given (and 
partially unknown) environment. The centrepiece of such 
an agent, which is also equipped with sensors (to receive 
perceptual input, percepts, from the environment) and 
actuators (enabling it to act on and change the envi¬ 
ronment) is a specific type of memory (ECM) [16]. In 
abstract terms, the memory is represented by a space of 
clips that can represent percepts, actions, and combina¬ 
tions thereof. After receiving sensory input, the PS agent 
initiates a random walk within the clip space to find an 
action. Under a given reward scheme, the agent’s choices 
have consequences that modify and update its memory, it 
learns. Throughout the learning process, the agent must 
hence be able to adjust its deliberation according to its 
experience, which entails updating the stochastic matrix 
P = (pij) of transition probabilities of the random walk 
in the space of memory clips. Recently, the PS model 
has been shown to perform competitively in typical arti¬ 
ficial intelligence benchmark tasks [53], and, further, that 
its memory structure provides a dynamic framework for 
generalization [54]. 

A particular variant of the PS that we shall focus on 
now, is reflecting projective simulation (RPS) [15]. In the 
RPS framework, the random walk in the memory space 
is continued until the underlying Markov chain P = {pij ) 
is (nearly) mixed, and actions are then sampled accord¬ 
ing to the stationary distribution tt, where Ptt = tt. The 
quantized version of the RPS yields a quadratic speed¬ 
up with respect to its classical counterpart, both in the 
number of calls to P needed to mix the chain, and in 
the number of samples until an action is obtained. The 
Szegedy-type quantum random walk involved in this pro¬ 
cedure requires several levels of coherent controlization. 
At the lowest level of this nested scheme of adding con¬ 
trol, one encounters a set of unitaries that en¬ 

code the n X n stochastic matrix P of an n-clip network. 
That is, the first column of Uj has real and positive en¬ 
tries {i = 1,2,... ,n), and each such probability uni¬ 
tary Uj may thus be parameterized by (n — 1) real an¬ 
gles Oj^i ,..., Qj,n-\. After each step of the learning pro¬ 
cess, the matrix P is updated according to the rewards 
that may have been incurred, requiring also updates of 
the Uj. In the spirit of adaptiveness of the agent’s design 
it is desirable that these updates can be carried out by di¬ 
rectly updating the angles 9j^k in an otherwise fixed hard¬ 
ware. A method [48] for realizing this requirement is a 
nested construction of coherent controlization for log 2 (n) 
qubits, where control is added to (n — 1) single-qubit Y- 
rotations U{9j k) for each Uj. 

For instance, for two-qubits, the unitary U{9i) is un¬ 
conditionally applied on the control qubit, followed by 
the applications of C/( 02 ) and U{93) conditioned on the 
control qubit being in the states |1)^ or |0)^, respec¬ 


tively, see Fig. 1 (b). For three qubits, a pair of two-qubit 
subroutines of the form just described replaces the two 
conditional single-qubit operations, see Fig. 3 (a), and so 
forth. This construction already entails the conditioning 
of single-qubit operations on all subspaces of the control 
qubits. We shall therefore consider the implementation 
of the probability unitaries Uj in superconducting qubits 
as a representative example that demonstrates the feasi¬ 
bility of coherent controlization. In order to proceed, we 
shall next give a brief overview of the properties of the 
superconducting qubit systems suitable for our purposes. 


C. Superconducting qubits coupled to microwave 
resonators 

The physical system that we consider in our proposal 
is an array of superconducting transmon qubits [31] cou¬ 
pled to a microwave resonator. For the reader unfamil¬ 
iar with the principal design of such a system, let us 
give an intuitive example. Consider a basic supercon¬ 
ducting LC-circuit, which may be thought of as the re¬ 
alization of a quantum mechanical harmonic oscillator, 
where charge and flux take the role of the canonically 
conjugate variables. Via the non-linearity of Josephson- 
junctions, an anharmonicity can be introduced into the 
system, which modifies the otherwise equal energy level 
spacing. This allows one to frequency-address transitions 
between two chosen levels (typically the two lowest-lying 
levels), thus forming a qubit. For the practical realiza¬ 
tion of such a macroscopic qubit, several options, such as 
charge and flux qubits, are available and we direct the 
reader to pedagogic reviews (see, e.g.. Ref. [29]) for more 
details on their differences. 

Here, we shall focus on the transmon qubit, introduced 
in Ref. [31], which is in its design similar to usual charge 
qubits, in the sense that two superconducting islands 
are connected via a Josephson junction with associated 
Josephson energy Ej. In some setups two Josephson 
junctions can be used to form a dc-SQUID. This leaves 
the possibility to modify the Josephson energy Ej of the 
junctions by threading the SQUID with an external mag¬ 
netic flux. Besides Ej, the energy levels of such charge 
qubits are determined by the charging energy Eq of the 
superconducting island. The departure of the transmon 
from other designs lies in the introduction of a large 
shunting capacitance parallel to the dc-SQUID, which 
drastically reduces Eq. The transmon qubit is hence op¬ 
erated in a regime where Ej Eq, which leads to an 
exponential decrease of the charge dispersion in Ej/Eq, 
while the anharmonicity is only diminished polynomially 
in this ratio. In other words, the qubit levels remain 
addressable by frequency selection, while their sensitiv¬ 
ity to environmentally induced charge noise is practically 
removed, which significantly improves qubit coherence 
times (> 100 (xs) [32, 33]. 

For our investigation, several such transmon qubits 
shall be considered to be capacitively coupled to a su- 
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perconducting resonator. To good approximation (see, 
e.g., the Appendix), the Hamiltonian for this dispersive 
interaction may be written as [55, 56] 

H/h= ujj.a)a + - y^ Xgjv ab\b^ 

i i 

( 2 ) 

i 

where and uiq^ are the (angular) frequencies of the 
resonator mode and the i-th qubit, respectively, with the 
corresponding ladder operators a,a\ b^,b\. For the re¬ 
mainder of this paper, higher excitation numbers are ig¬ 
nored for the qubit-modes. Due to the coupling of the 
qubits to the resonator, the latter also acquires an anhar- 
monicity, which causes the undesired self-Kerr effect [57] 
represented by Xrr- The constants Xrr^ Xqiqa the 
cross-Kerr coefficient Xgiv are hence not independent. For 
instance, when only one qubit is present, qi = q, the re¬ 
lation between these parameters is Xrr = X^rli'hCqq) in 
the dispersive limit [55]. Typical values for these pa¬ 
rameters that we will consider in the two-qubit case are 
Xg.r/(27r) « 1 —5 MHz, and Xgq/(27r) « 300 MHz (where 
we have assumed Xqiqi = Xq 2 q 2 = Xqq)i correspond¬ 
ing to cavity anharmonicities Xrr/(27r) roughly between 
0.7 - 17 kHz. 

The terms proportional to Xqir can be interpreted as 
conditional frequency shifts: Depending on the number 
of excited qubits, the resonator frequency ujr can be re¬ 
garded as shifted by the sum of the corresponding values 
XqiT- Conversely, the frequency of the i-th qubit can be 
considered as being shifted by the product of Xqir and the 
number of photons in the cavity. In the strong dispersive 
regime [45], the spectral lines of the cavity for different 
qubit states are well-resolved, and, vice versa, so are the 
qubit transition frequencies for different numbers of ex¬ 
citations of the resonator. Both of these points of view 
will play an important role in our protocol for coherent 
controlization that we will introduce next. 


III. COHERENT CONTROLIZATION EOR TWO 
SUPERCONDUCTING QUBITS 

In the protocol for coherent controlization that we pro¬ 
pose here, the phase space of the resonator mode serves 
as a bus between the qubits. It enables the conditioning 
of operations on particular subspaces of the qubit Hilbert 
space by separating the resonator states corresponding to 
different qubit subspaces in phase space. The mechanism 
for this separation is the free time evolution of coherent 
states with different frequencies. Recall that the reso¬ 
nance frequency of the cavity depends on the state of the 
qubits. Similar to the procedures used in [44, 45, 58], the 
following operations are employed for our protocol: 

• Unconditional displacements Da'. A very short 
pulse (a few ns), that is, sufficiently broad in fre¬ 
quency (of width ^ J^iXqir), so as not to distin¬ 


guish between the state-dependent frequencies of 
the resonator, displaces the cavity state indepen¬ 
dently of the state of the qubits. 

• Free time evolution U{t): During the free time 
evolution of the cavity mode, governed by the cor¬ 
responding parts of the Hamiltonian in (2), coher¬ 
ent states of the resonator corresponding to differ¬ 
ent qubit states rotate in phase space at different 
speeds. Appropriate waiting periods can hence be 
used to separate or recombine different coherent 
states. 

• Conditional qubit operations: When the (av¬ 
erage) photon number in the resonator is h, the 
(mean) qubit frequencies uig. are shifted to uiq. — 
Xqtrn with a spread XqirV^- Addressing the qubits 
with signals sufficiently narrow in frequency around 
LUq- therefore conditions the single-qubit operations 
on the cavity vacuum state JO)^. Unconditional 
qubit operations can be realized by appropriately 
broad (or multi-frequency) pulses. 

With these operations available we shall now specialize 
to the case of two qubits. 

A. Ideal two-qubit protocol 

Let us first analyze an idealized situation where the 
cavity Kerr effect due to Hx/fi = — and any 

loss due to interactions with the environment can be dis¬ 
regarded. Once we have set up our protocol for two 
qubits in this idealized setting, we shall consider the ro¬ 
bustness of the protocol under the influence of the men¬ 
tioned harmful effects. Irrespective of this restriction, the 
cross-Kerr interaction term Hi/h = — '^iXqir o-b\b- 
(in the strong dispersive regime) splits the cavity reso¬ 
nance frequency into 4 well-resolved spectral lines Wooj 
woi, wio, and cou, corresponding to the two-qubit ba¬ 
sis states [00)^, JOI)^, JlO)^, and 111)^. By selecting 
the qubit-cavity cross-Kerr coefficients Xq^rj^ = Xq 2 r = 
Aw « 1 — 5 X 27r MHz, the spectral lines are equally 
spaced, i.e.. 

Woo ~ ^01 “b Aw = wio -b 2 Aw = wii -t- 3 Aw = Wj., ( 3 ) 

which can be abbreviated to Wmn = uJr — (2m -I- n)Aw. 
With the resonator prepared in the vacuum state JO)^, 
we now construct a protocol that implements coherent 
controlization for a set of three single-qubit unitaries 
{C/(0i)}i=i_2,3 to realize the probability unitary 

Cb = (i0)(01,^ ® U{es)q, +11)( 11^^ ® U{e2)q,) 

X (f/(0i),2 , (4) 

represented by the circuit-diagram in Fig. I (b). The 
protocol, whose steps are detailed in Fig. 1 (a) and which 
are illustrated in Fig. I (b) and (c), can be decomposed 
entirely into the operations described above. 
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(a) Ideal two-qubit protocol 

(i) The first (unconditional) single-qubit unitary U{6i) 

is applied to qubit 2, mapping the initially (pure) 
state \ |0)^ to (coo | 00)^ -I- coi | 01)^ -I- cio 1 10)^ + 

cii 111 )g) |0)^, where the coefficients Cmn = Cmn{6i) 
depend on the initial state and U{9i). 

(ii) An unconditional displacement Da maps | 0)^ to \a.)^. 

(iii) In the frame rotating with ojoo = ^r, the coherent state 

components for the qubit states rotate with (angu¬ 
lar) frequencies 0, Atu, 2Aa;, and SAcu for | 00)^, | 01)^, 
I 10)^, and I ll)q, respectively. That is, the time evo¬ 
lution Ui = jg governed by the cross-Kerr 

interaction term Hi. A waiting period of At = -nj/S.u) 
hence results in the state (coo | 00)^ -|-cio 1 10 )g) | a),, + 
(coi|01),-bcii|ll)J| —a)^, separating the subspaces 
for I 0)„ and 11)„ . 

(iv) Another unconditional displacement Da shifts the co¬ 
herent state components for the subspace | 1)^^ lo the 
vacuum, |-a)^ i-s- |0)^, while |a)^ i->- |2a)^. 


(v) During another waiting period of duration At the sec¬ 
ond single-qubit unitary 11 ( 62 ) is applied conditionally 
on the cavity being in the vacuum. Meanwhile, the 
coherent state components |00 )^|q)^ and |10 )^|q)^ 
complete 0 and 1 revolution, respectively, returning to 
the initial position of step (v). 

(vi) An unconditional displacement by —2a exchanges the 
roles of the subspaces |0)^^ and 11)^^' 

(vii) During the penultimate waiting period of At, which ro- 
tates(|01)^ + |ll)J|-2a)^to(|01)^ + |ll)J|2a),, 
the third single-qubit operation 1 /( 63 ) is applied con¬ 
ditionally on the cavity vacuum state. 

(viii) An unconditional displacement D-a shifts |2a)^ to 
|a)^, and |0)^ to |-a)^. 

(ix) The last period of free time evolution by At rotates all 
phase space components to |—a)^, disentangling the 
cavity from the qubits. 

(x) The final unconditional displacement Da returns the 
cavity to the vacuum state. 


(b) 


(c) 




<x> o 100), • I01>, • lio>, • 111), 


FIG. 1. Two-qubit protocol for coherent controlization: The steps (i)-(x), listed in (a) realize the probability unitaries 
Uj (positive, real entries in the first column), represented by the circuit diagram in (b). The filled dots on the controlled 
operations indicate that the unitaries on the target are conditioned on the control qubit state 11) , while the hollow dots 
“o” denote conditioning on the control qubit state |0)^^, by way of coherent controlization of the three single qubit unitaries 
U(6i), 1/(62), and 1/(63). In (c), an extended circuit diagram is shown, where the upper and lower fines indicate qubits q2 
and q\, respectively, while the middle line labelled C represents the cavity mode. The white rectangles represent unconditional 
displacements, and the control lines between the cavity and the single-qubit operations indicate operations conditioned on the 
resonator being in the ground state. The waiting time between the displacements is At = tt/Aoj. (d) shows the phase space 
representation of the protocol for the coherent state components corresponding to different qubit states. 


Given that the initial (unconditional) qubit rotation 
U(9i) can be performed very rapidly (« 10 ns), the 
overall gate time T 2 QB is approximately T 2 QB ~ 4At = 
47r/Aw = 400 ns to 2 fxs [for Xq-irl^T^) between 5 MHz 
and I MHz], compared to typical [44] cavity coherence 
times Tr ~ 100 /iS, and qubit relaxation- and dephasing 
times of around 20 — 100 /xs. To maintain a high fidelity, 
especially when the protocol is extended to several qubits 
(as we shall do in the following), it is desirable to de¬ 
crease the overall gate time, and hence to increase the 
qubit-cavity cross-Kerr coefficients Xqiv- However, this 
can only be done at the expense of an increased cavity 
self-Kerr term proportional to Xrr-, which grows quadrat- 
ically with the Xqiv- In other words, there is a trade-off 
between decoherence and imperfections of the gates due 


to distortions of the coherent-state components (through 
the Kerr effect). Before we extend our protocol to larger 
numbers of qubits, we will therefore study the robustness 
of the two-qubit protocol under these effects. 


B. Influence of the cavity self-Kerr effect 

With increasing coupling between the qubits and the 
resonator, the influence of the anharmonicities of the 
transmon qubits on the resonator mode becomes ever 
stronger. The ensuing cavity Kerr effect distorts (and 
displaces) the shape of the phase space distributions [57] . 
To estimate the impact on the performance of our pro¬ 
tocol, we shall more closely inspect the three stages at 
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which the overlap of the phase space distributions with 
the target states determines the success of the protocol. 
First, let us denote the overall state after step k as 

■ (5) 

m,n—0yl 

During steps (iv) and (vi), the coherent state compo¬ 
nents IV'oT’V- \'<Pii)r and IV’ooOr. corresponding 

to the subspaces and |0) , respectively, are re¬ 

quired to significantly overlap with the vacuum to achieve 
the conditioning of 17(02) and U{9z)- Finally, after being 
disentangled from the qubits in step (ix), all resonator 
states (l,j = 0,1) should ideally return to the 

vacuum in step (x). The resonator states at these vari¬ 
ous stages can be obtained by applying the displacements 
and waiting periods as described in the Fig. 1 (a), where 
the time evolution is now determined by Hi + Hk, with 
HK/k,= i.e., 

li’mDr = DaUi.K{M)Da\0)^\mn)^ , (6a) 

IV'miOr = D_ 2 aUl^K{At) | V"));,) )r- \^^)q ^ (6b) 

IV'mUr \^n)q = D^Ulk{A t)D_^Ui^K(At) | \mn\ , 

(6c) 

where 17/,^(t) = Ui{t)UK{t) = and 

At = tt/Auj. Note that the operations Uj can be elim¬ 
inated from (6), since [Hi,Hk] = 0 and the displace¬ 
ments Da = exp(Q!a^ — a*a) satisfy exp(i(j)a^Da = 
exp(j^a^a). Acting on a state with fixed m,n, 
one therefore has 

Ui{At)Da = 

(7) 

where the rotating frame frequencies are = ujmn — 
LOr = —{2m + n)Auj. The [//(At) can hence be easily 
commuted with the displacements Da (at most changing 
the signs of the displacement parameters for (m, n) = 
(0,1), (1,1)). In the parameter regime that we have cho¬ 
sen, one further has At/2 « 6.5 x 10“^ —3 x 10“^ 1, 

where the couplings of the two qubits to the cavity con¬ 
tribute to the cavity Kerr term approximately (see the 
Appendix) as Xrr « {xI,t + xl^r)/i^Xqq)■ We may 
therefore expand t7i^(At) into a power series for small 
c = i.e., 

2 

UK{At) = = 1 +te(a^a)" - ^-{a}a)^ + 0{e^), 

( 8 ) 

where 0{x) is a quantity such that 0{x)/x is bounded 
in the limit a: —>■ 0. Since D_a = D/, = D~^, we can 
then use the simple identities D-aO-^Da = + a* and 

D-adDa = a + a. With some tedious but straightfor¬ 
ward algebra we arrive at the estimates for the desired 
overlap after step (iv), i.e., 

= 1 - + 6n^ + n) + 0{e^) , (9) 


where n = \a\’^ and we have used the shorthand 

= \q{mn\^{0\^li'^^)^\mn)^\^ . ( 10 ) 

Similarly, we find the overlaps after steps (vi), 

= 1 - e2(409fi3 i5g^2 9fj) +c>(e^), 

( 11 ) 

and for step (x) for all to, n = 0 , 1 : 

pp^ = 1 - e2(256fi3 -h 1367/2 -h An) + ©(e^). (12) 

The explanation for the revival of the fidelity from 
step (vi) to (x) lies in the leading order Kerr effect. As 
discussed in [44, App. B], to linear order in e, the Kerr ef¬ 
fect does not distort the shape of the coherent states, but 
adds an amplitude-dependent rotation and global phase 
factor to each coherent state, that is 

g*.(ata)= 1^)^ = ^ (jg) 

Although the phases appear as relative phases be¬ 

tween different qubit states, our protocol is designed in 
such a way that all qubit states acquire the same (and 
hence a global) phase: in steps (hi) and (ix) all compo¬ 
nents have the same amplitude jaj, and the phases picked 
up by 100 )^ and | 10 )^ in step (v) are compensated by 
those acquired by jOl)^ and 111 in step (vii). 

The additional rotation(s) from a to on the 

other hand, are only partially compensated in the proto¬ 
col of Fig. 1. The coherent state components [00)^ and 
110 )^ overshoot their target in steps (hi) and (v), i.e., 
the coherent states are rotated counter-clockwise with re¬ 
spect to their target on the horizontal axis. This means, 
the displacement D-a in step (vih) leaves them lagging 
behind, that is, rotated clockwise w.r.t. their target. 
This lag is partially compensated by the over-rotation 
in step (ix) in the sense that the mismatch in step (iii) 
is fully corrected, but the error incurred during step (v) 
remains. A similar argument applies for the components 
corresponding to | 10 )^ and [ 11 )^, resulting in an im¬ 
proved fidelity in (12) as compared to (11). Nonetheless, 
the relatively large influence of the Kerr effect in (9)-(ll) 
is problematic, and indeed does not justify terminating 
the power series after order e^. We shall hence modify 
our protocol in a similar fashion as discussed in [44] to 
correct for the leading order Kerr effect altogether. 

C. Corrected two-qubit protocol 

To completely compensate for the contribution of the 
linear order Kerr effect, the displacements after each 
period of free time evolution are adjusted by angles 
tp-y = e( 2 | 7 |^ -|- 1 ), where 7 is the maximal displacement 
of the coherent state components during the time evolu¬ 
tion. Alternatively, this may be seen as a suitable change 
of basis in the phase space, as illustrated in Fig. 2 (a). 
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FIG. 2. Corrected two-qubit protocol: (a) The ideal two-qubit protocol from Fig. 1 is modified to correct for the 
leading order cavity Kerr effect. After each period of free time evolution, the basis in the phase space is adjusted by an angle 
= e(2|ap + 1) [steps (hi) and (ix)] or ip 2 a = £(8|ap + 1) [steps (v) and (vii)], corresponding to the largest displacement in 
the previous step, (b) The gate fidelity (top three blue curves) as measured by J-mh (the overlap of the final cavity states with 
the vacuum) and the gate time T2 qb = 4At (in fis, bottom green) are shown as functions of Xq2r/{‘2'^) = Auu/{2tt) in MHz. 
The fidelities (not taking into account decoherence) are shown for values Xggli'^'’^) = 300 MHz, and average photon numbers 
h = 1.5 (solid, top), h = 2 (dashed, second from the top), and h = 3 (dotted, third from the top), (c) A simulation (using 
the QuTiP library [59]) of the Wigner function of the reduced cavity state in the corrected protocol steps (i)-(x) from (a) is 
shown for h = 3, Xg2r/{‘^'^) ~ 1-5 MHz, and Xgg/{‘^'^) = 300 MHz [corresponding to the third curve from the top in (b)] when 
dephasing and amplitnde damping for the qubits (with coherence times = r, = 100 /rs), and photon loss for the resonator 
{Tr = 100 fis) are taken into account. Note that even without loss the reduced cavity state n |cmrip j V’mnj does 
not feature interference fringes. In addition, the application of the conditional qubit unitaries can decrease the interference 
between the vacuum and excited state components in the reduced cavity state, partially restoring the rotational symmetry of 
the Wigner function. This can be seen, e.g., in the transition from step (vi) to (vii) in (c). The simulated protocol, for which 
the reduced qubit state amplitudes \q{^‘i'\ (Tirp) \mn}^\ = 0,1) are shown in (d), was run for the initial qubit state 

I 00)^, and with rotation angles 29\ = 82 = O3 = n. Including decoherence, the simulation with n = 3, Xg2rl{‘2Tf) = 1.5 MHz, 
and Xggli‘2'^) = 300 MHz yielded J-mh = 95% and the overlap of the final state of the qubits with the target state amounted 
to 93%. Further details on this, and additional simulations with other choices of parameters can be found in the Appendix. 


For the evaluation of the fidelities, we hence include 
the transformations Uq,{'y) = ex'p(—i(p.yd''a), that is 

l^mn)r = DaUq,{a)Ui^Ki^t)\0)^\mn)^ , (14a) 

l^mn)r 1 ^^)^ = (2a) C//,iy (At) | , 

(14b) 

l^mUr \'^^)q = 17^ (o) C//,iC (At)£)_a (2a) 

X Ui.xiAt) \mn)^ . (14c) 

Proceeding as before, we arrive at the corrected fidelities 
for the overlap with the vacuum after steps (iv) , (vi) , and 
(x), given by 

= 1 - 2n‘^e^ + 0{e^) , (15a) 

= 1 - > (15b) 

J-Xl = 1 - 72nh^ + Oie^) . (15c) 

The fidelity of the protocol as quantified by J^mh (not 


yet taking into account decoherence) hence also depends 
on the displacement via the average photon number n = 
jap. Crucially, these values have to be chosen such that 
the overlaps between the different coherent state compo¬ 
nents remain small in steps (v) and (vii) . For values fi = 
1.5 and 2 one finds |(0| 2a)p = exp(—4|ap) « 2.5 x 10“^ 
and 3.3 x 10“^, respectively. For Xqq/{‘^'^) = 300 MHz 
and n = 1.5, one may reach fidelities J^mh between 99% 
and 94% for qubit-cavity cross-Kerr coefficients ranging 
from Xq 2 r-li‘^T^) = 1 MHz to 3 MHz, see Fig. 2 (b). In 
the Appendix we further include decoherence effects — 
dephasing and amplitude damping of the qubits with co¬ 
herence times T 0 = 30 — 100 /xs and r, = 20 — 100 /xs, 
respectively, as well as photon loss in the cavity with 
coherence times = 100 /xs, but we consider the single¬ 
qubit operations to be perfect. Employing simulations 
coded in PYTHON using the QuTiP library [59] we find 
that fidelities of 95% are reasonably achievable, e.g., for 
h = 3, Xq2r/(27r) = 1.5 MHz, Xqq!^'^) — 300 MHz and 





















































coherence times of 100 /iS for a range of initial states and 
angles 9i [i = 1, 2, 3), e.g., as shown in Fig. 2 (c,d). 


D. Other corrections 

Before we finally turn to the extension to three (and 
more) qubits, let us remark on additional possible sources 
of decreases in fidelity. In a similar way, in which the 
cavity inherits the Kerr term from the anharmonicities 
of the qubits, the qubits are also coupled via the cavity, 
even if direct capacitive coupling can be avoided by ar¬ 
ranging the qubits to be spatially well-separated. The 
corresponding term of the form Xqiq 2 ^\^i ^\^2 leads to 
additional phases that are acquired by the components 
IV'n )r 1 11)5 during the time evolution. This effect is at 
most of a size comparable to the cavity Kerr effect, which 
is required to be small to achieve high-fidelities. 

However, when the effect becomes non-negligible, it 
can be corrected by applying an appropriate phase gate 
along with 17(02) during step (v). In a regime where this 
becomes necessary, additional measures are further re¬ 
quired to compensate the cavity Kerr effect beyond the 
linear order corrections that we have considered so far. 
If required, this can be achieved by a scheme, recently 
proposed in Ref. [60], that relies on additional ancilla 
qubits. In the strong dispersive regime, the frequency 
shift of the ancilla qubit due to different photon numbers 
in the cavity is used for photon-number selective phase 
gates that compensate the phases acquired by the differ¬ 
ent Fock state components of the coherent states. When 
such corrections are used to eliminate the cavity self-Kerr 
effect, the fidelity iFmn of our protocol with parameters 
as in Fig. 2 can reach 99%. 

Another potential source of errors lies in imperfections 
in the circuit fabrication that may cause a small devia¬ 
tion e from the desired ratio Xqir = of the qubit- 
cavity cross-Kerr coefficients, that is, one may encounter 
a situation where Xqir = + for some small (real) 

I (51 ^ 1. This causes an additional rotation for the coher¬ 
ent state components corresponding to |01)g and |ll)g 
in steps (iii), (vii) and (ix), but not for the components 
corresponding to | 00)^ and 110)^. For steps (iii) and (vii) 
this effect can be compensated by modified displacements 
in steps (iv), (vi) and (viii). However, to recombine all 
coherent state components in step (ix), an “echo”-type 
operation is required. That is, step (ix) is amended in the 
following way: after a waiting time t' = ^^(fqif) < At 

an unconditional 7r-pulse [44] is applied to the second 
qubit, which exchanges the subspaces {|00)^,|10)^} and 
{|01)g,|ll)^}. After another waiting period of duration 

t" = — -1' a final 7r-pulse to the second qubit restores 

the original qubit state. Step (x) can then be executed 
with an appropriately modified displacement to complete 
the protocol. With the possibility for these corrections 
in mind, we now turn to the three-qubit protocol. 


IV. COHERENT CONTROLIZATION FOR 
THREE QUBITS & SCALABILITY 

Let us now discuss the extension of our protocol to 3 
qubits and beyond. We will consider the ideal-three qubit 
protocol for the realization of a probability unitary Uj as 
a proof-of-principle example from which the construction 
of the protocol for any number of qubits may be inferred. 
To extend the previous protocol, first note that the cross- 
Kerr interaction Hi/h = — '^iXqiv 0 ‘b\b^ again splits 
the resonator frequency into spectral lines for all the 
qubit states. The sidebands are equally spaced by Aw by 
selecting Xqir/'^ = Xqir!^ = Xq^r = Aw. The ideal pro¬ 
tocol, detailed in Fig. 3, then relies on the same type of 
operations as before, unconditional displacements (now 
with amplitudes up to 4a), and waiting periods (now of 
durations At/2, At, and 2At), and qubit operations con¬ 
ditioned on the cavity being in the vacuum state. 

Increasing the number of qubits hence brings about 
exponential scaling of several quantities of interest — as 
expected when exponentially increasing the dimension of 
the state space that we wish to explore. To separate the 
coherent state components corresponding to all different 
subspaces, larger displacements and finer graining of the 
rotation angles in phase space are required. That is, in 
addition to rotations by half-periods, for three qubits also 
quarter periods are necessary, and for n-qubits, rotations 
by 27r/2"'“^ are needed. In addition, waiting periods 2At 
corresponding to the full period become obligatory when 
applying the conditional qubit rotations. Starting from 
an ideal n qubit protocol (for n > 3), an additional qubit 
may be added by inserting 2"“^ sequences of operations 
between each pair of conditioned unitaries on the qubit 
previously labelled qi (e.g., C/( 04 ) and U{9^) for n = 3). 
Each sequence consists of (at most) 4 displacements and 
two conditioned unitaries on the new qubit and the wait¬ 
ing periods after the 4 new displacements are of duration 
At/2"“^, 2At, 2At and 2At(I — 2“”), respectively. We 
hence find that, in an ideal n-qubit protocol, the number 
of displacements is N£,(n > 3) = 2"+^ — 4, the number 
of qubit unitaries is Nu{n) = 2"’ — I, and the overall 
duration is Tin > 3) = (3 x 2" — II)At. In a nonideal 
protocol, correction operations have to be included as 
well. For instance, for each conditioned qubit unitary 2 
echo pulses (as discussed in Section HID) may need to 
be added, increasing the total number of qubit opera¬ 
tions to 3 X 2”“^, but leaving the gate time unchanged. 
The exponential increase in dimension hence carries over 
to the scaling of the number of conditional qubit oper¬ 
ations, and to the overall gate time. Since increases in 
gate times, displacements and qubit-cavity cross-Kerr co¬ 
efficients XqiT all lead to increased disturbance due to the 
cavity Kerr effect, its compensation using photon-number 
selective phase gates [60] becomes ubiquitous, despite the 
possibility to correct for the linear order Kerr effect in 
Eq. (13) even in the presence of three individual phase 
space components with different displacements. 

Nonetheless, it should be mentioned that such chal- 
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FIG. 3. Three-qubit protocol for coherent controlization: The circuit shown in (a), where the filled dots on the 
controlled operations indicate that the unitaries on the target are conditioned on the control qubit state 11)^^, while the hollow 
dots “o” denote conditioning on the control qubit state lO)^^, represents a three-qubit probability unitary Uj. (b) shows an 
extended circuit diagram that realizes the circuit in (a) by adding control to the single-qubit y-rotations U{ 9 i). The uppermost, 
and the two lowest horizontal lines indicate qubits 3, 2 and 1, respectively, while the second line from the top represents the 
cavity mode. The white, diamond-shaped rectangles represent unconditional displacements of all cavity modes, and the control 
lines between the cavity and the single-qubit operations indicate operations conditioned on the resonator mode being in the 
ground state. The waiting time between the displacements varies between At/2 = 7r/(2Auj) and the full period 2At. (c) shows 
the phase space representation of the coherent state components (*,/, fc = 0,1) corresponding to the qubit states | ijk) 

for the steps I = (i), , (xxiv) of the ideal three-qubit protocol. 


lenges for the scaling of quantum computational architec¬ 
tures are expected for every physical platform, and are 
hence not specific to our proposal. For our three-qubit 
protocol, simulations (using the QuTiP library [59]) that 
take into account partial correction of the linear Kerr ef¬ 
fect and decoherence yield fidelities for the cavity state 
of 80% (64% for the qubits) for the system parameters 
Xq3r/(27r) = 0.3 MHz, Xqqli‘^'^) = 300 MHz and h = 1, 
see the Appendix. However, when corrections as dis¬ 
cussed in [60] are included, the same parameters yield 
fidelities of up to 95% for the cavity, and 75% for the 
qubits. 


V. DISCUSSION 

We have introduced a protocol for coherent controliza¬ 
tion, that is, adding control to (a set of) unspecified or 
unknown unitaries, using superconducting qubits cou¬ 
pled to a microwave resonator. This task is of interest for 
the flexible realization of quantum computational archi¬ 
tectures, but also of great importance for the adaptive¬ 
ness of quantum-enhanced learning agents [48] . We have 
selected an example from the latter context, the reflec¬ 
tive projective simulation model [15] for artificial intelli¬ 
gence, where coherent controlization is already useful at 
the lowest level of the deliberation algorithm to construct 
coherent encodings of Markov chains — the probability 
unitaries. 
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We have given explicit protocols for the realization 
of these unitary operations using transmon qubits [44], 
Within the strong dispersive regime, we exploit the cou¬ 
pling between the qubits and the resonator. The cav¬ 
ity mode here acts as a bus between the qubits, playing 
the role of the additional degree of freedom necessary for 
adding control to unknown unitaries, similar, e.g., to the 
vibrational modes of trapped ions [51]. We have provided 
a detailed discussion of the role of the cavity Kerr effect, 
the strongest source of disturbance, in our protocol, in¬ 
cluding corrections for the linear order effect. Based on 
these considerations, and bolstered by numerical simu¬ 
lation including photon loss, as well as amplitude- and 
phase-damping for the qubits (featured in the Appendix) , 
we conclude that a possible experimental realization of 
our protocol with two qubits may achieve high fidelities 
(up to 95% for the cavity and 93% for the qubits) for 
reasonable ranges of the system parameters. We hence 


consider our proposal for two qubits to be readily imple- 
mentable using current superconducting technology. For 
three (or more) qubits, an implementation is still possi¬ 
ble, although the significant drop in fidelity (< 80% and 
< 64% for the cavity and qubits respectively) suggests 
that upscaled versions of our protocol may require addi¬ 
tional corrections of the Kerr effect [60] , which can bring 
the fidelities up to 95% (cavity) and 75% (qubits). 
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APPENDIX: 

SUPPLEMENTARY INFORMATION 

The Appendix is structured as follows. In Appendix A 
we provide a derivation of the Hamiltonian of the trans- 
mon qubits coupled to a resonator. The aim of this sec¬ 
tion is to give a transparent account of all approxima¬ 
tions made in the system we consider and to give an 
introduction to this kind of superconducting qubits that 
is accessible for non-specialists. In Appendix B we then 
present the details of the simulations that were conducted 
to assess the robustness of our protocol in the presence 
of the Kerr effect, dephasing and amplitude damping of 
the qubits, and photon loss in the resonator. 


Appendix A: Transmon qubits coupled to a 
microwave resonator 


A.I. Prom Josephson junctions to transmons 


A superconducting LC-circuit may be thought of as a 
harmonic oscillator, where the position and conjugate 
momentum variables are the flux $ through the in¬ 
ductor and the charge Q on the capacitor plates, re¬ 
spectively. Indeed, these quantities have to be treated 
as operators satisfying the canonical commutation rela¬ 
tion [$ , Q] = ih, see, e.g., Ref. [29]. In other words, the 
system Hamiltonian H = Q'^l{2C) + l{2L) may be 

written as 


H = huj{q^q J- 1/2) 

in terms of the ladder operators defined by 



(A.l) 


(A.2a) 


(A.2b) 


with ct) = I/'/LC, Z = sjLjC, and = 1- The 

introduction of a Josephson junction (a thin insulating 
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barrier separating two pieces of superconducting mate¬ 
rial) creates a non-linearity in the system that allows 
identifying two energy levels as the qubit levels. The 
Hamiltonian is then modified to 


H = 


2 Cs 


$2 

— — COS 



(A.3) 


where g is the qubit-cavity coupling, ujn and ujq are the 
frequencies of the isolated qubit and cavity, respectively, 
and the operators c and cl are the resonator ladder op¬ 
erators satisfying [c,cl] = 1. The terms in Hqr that are 
quadratic in the mode operators can be diagonalized by a 
Bogoliubov transformation. That is, one introduces the 
dressed operators a and b, such that 


where the junction inductance Lj is typically expressed 
via the Josephson energy Ej, that is, Lj = ^q/Ej, with 
the magnetic flux quantum $o = hf{2e), and ^/d’o is 
the phase difference across the junction. The original 
capacitance and the Josephson junction now form a su¬ 
perconducting island — a Cooper-pair box (CPB) — 
with capacitance Cs = C + Cj^ which can be written 
in terms of the charging energy Eq = I{2Cy,)- The 

charge in the CPB depends on n, the number of trans¬ 
ferred Cooper pairs (with charge 2e), and the effective 
offset charge 2eno, he., Q/(2e) = n — Uo- With the re¬ 
placement Z —>■ Zj = \JEj ICy:, one may dehne raising 
and lowering operators in full analogy to (A.2). To re¬ 
move the strong dependence of the energy-level splitting 
on the (environmentally induced) offset charge, an addi¬ 
tional shunting capacitor with large capacitance Cb can 
be inserted into the circuit in parallel with the Josephson 
junction to create the transmon [31]. This strongly in¬ 
creases Cy ^ Cy = Cb + C + Cj, such that Ec/Ej 1. 
Inserting the ladder operators q and ql (with Z —)■ Zj) 
into Eq. (A.3) and expanding in powers of Ec/Ej we 
arrive at 


c = cos 0 a + sin 9 b , 
q = — sin 9 a + cos 9 b , (A.7) 

with tan 20 = 2glA, where A = loq — ujr is the qubit- 
cavity detuning. The excitations of the dressed modes 
can no longer be uniquely attributed to just the res¬ 
onator or just the qubit. However, in the strong disper¬ 
sive regime, where g/A <?; 1, the mixing angle is small, 
9 Ki g/A and excitations created by (ql) are associ¬ 
ated “mostly” with the cavity (qubit). The Hamiltonian 
transforms to 


Hqr/h = ujrO^a + ujqb^b — 2j3 siv?9 cos^9 a'ab^b 

— — cos^9 b^b^bb — — sin^da^a^aa 
2 2 

-I- /3sin0cos0 sin^0(a^a^ab 4- a^aab^) 

+ cos^9 (^a^b^bb+ ab^b^b^ 

+ sin 0 cos 0 (^a'a^bb+ aab^b^')~\ , (A.8) 


H = ^ 2 EjEc -Ej-^ + (^SEjEc - Ec) q^q 

- ^ + q^qqq) - § {qqqq + q^^ q^ q^) 

-^(gg + gV + ^Vqq). (a.4) 

Within the rotating wave approximation, the Hamilto¬ 
nian of Eq. (A.4) can be truncated to the effective trans¬ 
mon qubit Hamiltonian 

Hq = hwQ q^q - ^q^q’^qq , (A.5) 

where wg = {\/8EjEc — Ec) /h is the transmon qubit 
frequency and /? = Ec/h is the anharmonicity. 


A.II. Transmons coupled to a resonator 

Next, we are interested in describing a transmon qubit 
that is capacitively coupled to a microwave resonator. 
Within the rotating wave approximation the joint system 
is described by the Hamiltonian 


where the dressed mode frequencies are given by 

uJr = ujrcos'^9 + u)Qsiv?9 — gsin(20), (A.9a) 

Aq = a;ijsin^0 -|- ujq cos^9 + gsin(20). (A.9b) 


The terms in the square bracket can be seen to oscillate 
rapidly, and we may therefore remove these terms in an¬ 
other rotating wave approximation. With the notation 
Xqr = (^/2) sin^(20), Xqq = P cos^0, and 

Xrr= = Psin^9, (A.IO) 

^Xqq 

we then arrive at the Hamiltonian of a single transmon 
qubit coupled to a resonator in the dispersive limit 


E[qj./h = ujrOpa -f ujqb^b 


Xrr 


aP aa 


= UJrOp a -I- UJqb^b 



^^b^ b^ bb — XqrCi^cbb^b 


^(b'^b)^ - XqTa^ab''b 

(A.ll) 


Hqr/h = wrc'c + wgq'l'q 4 - g {o'q 4- cq^) - ^q^q"'qq , 

(A.6) 


where uJr = + Xrr/2 and Uq = ujq + Xqq/‘2- The dressed 
qubit anharmonicity Xqq: the qubit-cavity cross-Kerr co¬ 
efficient Xgr, and the cavity self-Kerr coefficient Xrr can 
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be expressed directly via the coupling strength g and the 
detuning A as 


Xqq — 


Xqr — 


Xrr — 


^ f 1 + 1^1 ] 

4 ’ 


2 A 2 + V 


1^1 

4 I ^A2 + 4g2 


(A.12a) 

(A.12b) 

(A.12c) 


and the possibility for correcting it entirely using photon- 
number selective gates [60]. In addition, all simulations 
assume the presence of amplitude and phase damping 
for the qubits, and photon loss in the resonator. The 
unconditional displacements can safely be assumed to be 
perfect. We further assume the single-qubit operations 
to be perfect, given that the pulses are slow enough to 
address only the zero-photon subspace, but short enough 
to fit within the At time intervals of our protocol. 


B.I. Setup for simulations 


A.III. Two transmons coupled to a cavity 

When two transmon qubits are coupled to the cav¬ 
ity, we may write an analogue expression to Eq. (A. 6 ). 
That is, in the rotating wave approximation we have the 
Hamiltonian 


Hggr/fi = UJR c^c + UJQ^ q\qi + qlq2 
+ 9i (c^9i + cql) -b 52 (c^92 + 

- yglqWigi - y929^292 , (A.13) 

where we have neglected any direct coupling of the qubits 
to each other. In the two-qubit case (and beyond), an 
analytical diagonalization of the harmonic part (terms 
quadratic in the mode operators) becomes infeasible. 
However, using numerical methods and following simi¬ 
lar arguments as presented in the previous Sec. A.H one 
arrives at the effective Hamiltonian 

Hgqr/h= OJra'^a + '^uiq.b\b- - y^Xqira^a'b\b- 

i i 

i 

Assuming that the transmons can be fabricated such that 
Xqiqt = Xqq qubits, we have numerically checked 

that the cavity self-Kerr coefficient Xrr can be approxi¬ 
mated as 


Xrr 


x\ 


qir 


'~q2r 


4Xa 


(A.15) 


Appendix B: Simulations of coherent controlization 
using transmon qubits 

In this section we present the simulations of our proto¬ 
col for coherent controlization with two and three qubits. 
In the simulations, which were coded in PYTHON using 
the QuTiP library [-59], we discuss the influence of the 
cavity Kerr effect, including the (partial) correction of the 
linear Kerr effect by way of reference frame adjustments. 


For the simulation we truncate the Hilbert space of the 
resonator to be spanned by the Fock states |n)^ of pho¬ 
ton numbers n = 0,1,..., 100. This is a good approxi¬ 
mation since the maximal average photon number of the 
coherent states in our protocol is fimax = |2ap = 14 
{nmax = |4ap = 32) photons for 2 (3) qubits and the 
overlap with photon numbers larger than 100 is hence 
negligible. For the Hilbert space of the transmon qubits 
we will each only consider the lowest two eigenstates, 
i.e., the qubit levels. The effective Hamiltonian in the 
frame rotating with ujr and LUg for the resonator and qubit 
Hilbert spaces, respectively, is given by 

Hi^K/b= [Hi + Hk)/^ 
i 

= - '^Xq^ra)aaf(T~ - ^{a)af , (B.I) 

i 

where = crJ ± iaf are the raising/lowering operators 
of the i-th qubit. In addition to the free time evolution 
we will include corrections for the linear order of the Kerr 
effect by including rotations U^{x) = exp{—i(pja^a) after 
each period of time evolution, where 7 is the maximal 
displacement of the different coherent state components 
in the preceding step of the protocol. The conditional 
single-qubit operations will be represented by 

\0){0l(^Uq{9,) + ( 1 ,- |0)(0|^) 01,, (B.2) 

where the single-qubit Y-rotations on the second (and 
third) of the two (three) qubits are realized by a time- 
dependent drive U{6j) = exp(—The dura¬ 
tions T of these drives are taken to be the durations of 
the waiting periods in the corresponding protocol steps 
(T = At for two qubits). On top of the unitary time 
evolution, unconditional displacements, and conditional 
qubit operations, all of which will be assumed to be per¬ 
fect, we will consider decoherence in the system. In par¬ 
ticular, we assume that the dynamics of the overall state 
p of the joint cavity-qubit system during the waiting pe¬ 
riods of the protocol is governed by the master equation 
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dp 
dt 

+ ^ p - ) 

^ i 

+ (^-3^)S(-f^-|-rf- (B-3) 

As the output of interest of the simulations we consider 
the hdelities Jv for the resonator, i.e., the squared over¬ 
lap with the vacuum at the final step of the protocol, 
given by = r(0| (Tr,jp) |0)^, and Tq for the qubits, 
i.e., the squared overlap with the target state \ 4’)q, given 
by = q{'tp\ (Tr^p) \'4’)q- We evaluate these two fideli¬ 
ties separately because the success of the protocol is ul¬ 
timately only determined by achieving a large overlap of 
the reduced qubits state with the desired target state, 
irrespective of the final cavity state. On the other hand, 
the cavity fidelity gives a more detailed overview of the 
errors incurred by the Kerr effect in the resonator. There¬ 
fore, both and JFq are quantities of interest for our 
protocol. In addition to the final fidelities, we include 
plots of the Wigner function W{x,y), where 

OO 

2 I (Top) |2:+ 2)r > 

— OO 

of the reduced resonator state pr = Tr^p, where |a;) are 
the eigenstates of the quadrature operator $ with eigen¬ 
value X, and plots of the absolute values of the density 
matrix elements Pq |m.n)^| of the reduced qubit state 
Pq = Tr^p throughout the steps of our protocol. 


i 1 

- [i^/,ic,p] -b — {2apa^ - a)ap - pa^a) 


B.II. Simulations for two qubits 

For two qubits, the qubit-cavity cross-Kerr coefficients 
Xqr = Xq^r = Xq\rl^ are chosen in a range from 
1.5 X 27r MHz to 3 x 27r MHz. According to the approxi¬ 
mation in Eq. (A.15), we further set Xrr = ^'}dqr!^^Xqq)i 
where Xqq is fixed to 300 x 27r MHz. The coherence 
time for the resonator is set to = 100 ps through¬ 
out. For the simulation results presented in Fig. B.l, we 
hence vary the qubit-cavity coupling, the displacements 
in the protocol, the decoherence times, the initial qubit 
state Pq = |C)(Clq and the angles 9j (j = 1,2,3) for the 
single-qubit operations. In addition to the results shown 
in Fig. B.l, we have simulated the influence of decoher¬ 
ence on the protocol when assuming that the cavity Kerr 
effect has been fully corrected using the methods of [60]. 
In that case, using the same parameters as for the simu¬ 
lation shown in Fig. 2 (c,d) of the main text, but setting 
Xrr = 0, we obtain = 99% and Tq = 96%. 

Finally, note that the parameters Xqtr, Xgtqi and Xrr 
are derived from the quantities Ec^, Ej., ujr, and the 
cavity-qubit couplings pi for i = 1,2. Achieving the 
exact required ratio Xqir = ' 2 'Xq 2 r hence relies on the 


fine tuning of these quantities. For instance, to ap¬ 
proximately achieve the parameters of the simulation 
in Fig. B.l (c), one may set Ej^ = 27.0 GHz, Ej^ = 
20.1 GHz, Ec, = Ec^ = 0.3 GHz, pi/(27r) = 101 MHz, 
P2/(27r) = 127 MHz, and ojR/{2Tr) = 9.16 GHz to ob¬ 
tain Xq 2 rli‘^'’^) = 1.499 MHz, Xqir/i‘^T^) = 2.982 MHz = 
1.989 Xq 2 r, X9i«i/(27r) = 296.998 MHz, Xq 2 q 2 li‘^'^) = 
298.46 MHz, and Xrrli^^T^) = 9.367 kHz (as compared 
to Xqir/i^Xqiqi) = 9-34 kHz). For the simulations, we 
assume that such deviations are corrected by an echo- 
type operation as explained in Sec. HID of the main text. 
Finally, note that all parameters used in the simulations 
are compatible with the dispersive approximation. 


B.III. Simulations for three qubits 

For three qubits, the qubit-cavity cross-Kerr coeffi¬ 
cients Xqr = Xqsr = ‘^Xq 2 r = ^Xq^r are choscn in a 
range from 0.2 x 27r MHz to 0.4 x 27r MHz. Accord¬ 
ing to the approximation in Eq. (A. 15), we further set 
Xrr = 21Xqr/(^Xqq)! where Xqq IS fixed to 300 X 27r MHz. 
The coherence times for the resonator and the qubit (de¬ 
phasing and amplitude damping) are set to 100 ps. For 
the simulation results presented in Table I, we hence vary 
the qubit-cavity coupling, the displacements in the proto¬ 
col, the initial qubit state pq = 10 (C Iq and the angles 9j 
{j = 1,2,... ,7) for the single-qubit operations. A selec¬ 
tion of these simulations Figs. B.3 and B.4 further illus¬ 
trate the reduced states of the resonator and the qubits 
throughout the protocol. 


# 

Xqr 

MHz] 

a 

\Oq 

01 

02 

03 

04 

05 

06 

07 

Tr 

Eq 

Si 

p 

CO 

X 

27r 

1 

\^u) 

0 

0 

0 

0 

0 

0 

0 

0.671 

0.296 

S2 

p 

CO 

X 

27r 

1 

jOlO) 

TT 

2 

TT 

0 

0 

TT 

0 

TT 

2 

0.505 

0.449 

S3 

0.4 

X 

27r 

1 

jOOO) 

0 

0 

0 

0 

0 

0 

0 

0.832 

1.000 

S4 

p 

CO 

X 

27r 

1 

jOOO) 

TT 

2 

0 

TT 

6 

0 

TT 

0 

TT 

2 

0.699 

0.481 

ss 

0.4 

X 

27r 

1 

jOOO) 

TT 

2 

0 

TT 

6 

0 

TT 

0 

TT 

2 

0.578 

0.457 

S6 

0.2 

X 

27r 

1 

jOOO) 

IT 

2 

0 

TT 

6 

0 

TT 

0 

TT 

2 

0.747 

0.444 

S7 

0.35 

X 

271 

1 

jOOO) 

IT 

2 

0 

TT 

6 

0 

TT 

0 

TT 

2 

0.646 

0.475 

S8 

p 

CO 

X 

271 

1 

jOOO) 

TT 

0 

0 

0 

0 

0 

TV 

0.697 

0.797 

S9 

p 

CO 

X 

271 

/ 1 
V 2 

jOOO) 

TT 

0 

0 

0 

0 

0 

TV 

0.674 

0.729 

SlO 

p 

CO 

X 

271 

fs 
V 2 

jOOO) 

TT 

0 

0 

0 

0 

0 

TV 

0.329 

0.796 

Sll 

p 

CO 

X 

271 

1 

jOOO) 

TT 

12 

0 

TT 

3 

0 

0 

TT 

2 

TV 

0.804 

0.639 

Sl2 

p 

CO 

X 

271 

1 

jOOO) 

TT 

12 

0 

TT 

3 

0 

0 

TT 

2 

TV 

0.953 

0.748 


TABLE I. Three-qubit simulations: The table shows the 
results of the numerical simulations for three qubits for vari¬ 
ous system parameters, rotation angles 9j, and initial states, 
where | ^„) is the uniform superposition overall three qubit 
computational basis states. The last simulation, Si 2 was exe¬ 
cuted for the same parameters as sn, but assuming that the 
cavity Kerr effect has been corrected (here: Xrr = 0) using 
the photon-number selective phase gates [60]. 
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# 

Xqr [MHz] 

a 

Tq [^S] 

l/rs] 

10, 

01 

O2 

^3 

Tr 


(c) 

3 X 27r 

y/2 

20 

30 

|00) 

TT 

2 

TT 

TT 

0.907 

0.916 

(d) 

1.5 X 271 


100 

100 

|00) 

TT 

2 

TT 

6 

TT 

3 

0.938 

0.953 

(e) 

1.5 X 277 


100 

100 

|00) 

0 

0 

0 

0.943 

1.000 

(f) 

1.5 X 277 


100 

100 

|00) 

TT 

12 

TT 

4 

TT 

0.955 

0.953 

(g) 

1.5 X 277 


100 

100 

|10) 

TT 

4 

TT 

TT 

6 

0.943 

0.959 


(a) 


(b) 


O 






FIG. B.l. Two-qubit simulations: (a) Ideal protocol, (b) Table of simulation results, (c)-(g) show the reduced resonator 
state Wigner function on the left-hand side, and the absolnte values of the density matrix elements \q{^^'V\ pq\m'n)^\ of the 
reduced qubit states on the right-hand side. 
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FIG. B.2. Three-qubit simulation si: (a) Ideal protocol, (b) The Wigner function of the reduced resonator state of the 
three-qubit protocol is shown for the parameters of the simulation specified in si of Table I. (c) shows the corresponding plots 
of the absolute values of the density matrix elements \q{iJLv\\ Pg \mnl}J {p, u, A, m,n,l = 0,1) of the reduced qubit states. 
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FIG. B.3. Three-qubit simulation sn: (a) Ideal protocol, (b) The Wigner function of the reduced resonator state of the 
three-qubit protocol is shown for the parameters of the simulation specified in sn of Table I. (c) shows the corresponding plots 
of the absolute values of the density matrix elements |g(/ri'A| pg \mnl}J {p, u, A, m,n,l = 0,1) of the reduced qubit states. 
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FIG. B.4. Three-qubit simulation S 12 : (a) Ideal protocol, (b) The Wigner function of the reduced resonator state of 
the three-qubit protocol is shown for the parameters of the simulation specified in S 12 of Table I, where the cavity Kerr effect 
is assumed to be fully corrected, (c) shows the corresponding plots of the absolute values of the density matrix elements 
|q(/rz 2 A| pq \mnl)^\ {p, v, A, m,n,l = 0,1) of the reduced qubit states. 







































































































